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ABSTRACT 



^ , We report sensitive ATCA radio continuum observations toward IRAS 15596—5301 

^ \ and 16272—4837, two luminous objects {C > 2 x IO^Lq) thought to represent massive 

star forming regions in early stages of evolution (due to previously undetected radio 
emission at the la level of 2 mJy per beam). Also reported are 1.2 millimeter continuum 
and a series of molecular line observations made with the SEST telescope. 

The radio continuum observations toward IRAS 15596—5301 reveal the presence of 
three distinct compact sources, with angular sizes of 2.7" to 8.8" (FWHM), all located 
within a region of 30" in diameter. Assuming that these are regions of ionized gas, we 
find that they have diameters of 0.06-0.2 pc, electron densities of 8 x 10^ — 2 x 10^ cm~^, 
and that they are excited by early B type stars. The 1.2-mm observations show that 
the dust emission arises from a region of 42" x 25" (FWHM) with a total flux of 5.8 Jy, 
implying a mass of 1.4 X 10^ Mq. The hue observations indicate that IRAS 15596-5301 
is associated with a molecular cloud with a FWHM angular size of 37" (~ 0.4 pc radius 
at the distance of 4.6 kpc), a molecular hydrogen density of ~ 4 x 10^ cm~^ and a 
rotational temperature of ~ 27 K. We suggest that the massive dense core associated 
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with IRAS 15596—5301 contains a cluster of B stars which are exciting compact H II 

regions that are in pressure equilibrium with the dense molecular surroundings. 

No radio continuum emission was detected from IRAS 16272—4837 up to a 3a limit 
of 0.2 mJy. However, the 1.2-mm observations show strong dust emission arising from 
a region of 41" x 25" (FWHM) with a total flux of 13.8 Jy, implying a mass of 2.0 x 10^ 
Mq. The line observations indicate the presence of an elongated molecular cloud with 
FWHM major and minor axes of 61" and 42"(0.50 x 0.35 pc in radius at the distance of 
3.4 kpc), a molecular hydrogen density of ~ 2 x 10^ cm^'^ and a rotational temperature 
of ~ 27 K. The high luminosity (2.4 x IO^Lq) and lack of radio emission from this 
massive core suggest that it hosts an embedded young massive protostar that is still 
undergoing an intense accretion phase. This scenario is supported by the observed 
characteristics of the line profiles and the presence of a bipolar outflow detected from 
observations of the SiO emission. We suggest that IRAS 16272—4837 is a bona-fide 
massive star forming region in a very early evolutionary stage, being the precursor of 
an ultra compact H II region. 

Subject headings: H II regions — ISM: individual (IRAS 15596-5301, IRAS 16272-4837) 
— stars: formation - stars: massive 

1. Introduction 

The earliest phase of high-mass star formation is possibly one of the least understood stage of 
evolution of massive stars. Massive stars (M> 8 Mq) are known to be formed in dense molecular 
cores, however the sequence of processes leading to their formation is not yet well established. In 
particular, the role of coalescence (Stabler et al. 2000) and accretion (Osorio, Lizano, & D'Alessio 
2000) processes in the assembling of a massive star is still under debate. The determination of the 
physical conditions of the gas during the formation and early evolution of a massive star is difficult 
because of their scarcity and rapid evolution. In addition, massive stars are usually born in clusters 
or groups hence their individual studies are usually afflicted by confusion, particularly because 
they are found located, on the average, at larger distances from the Sun than sites of low-mass star 
formation. Massive objects in early evolutionary stages, namely in the process of assembling the 
bulk of their final mass, should be identified by having high bolometric luminosities (> IO'^Lq), 
strong dust emission, and very weak or no detectable free-free emission at cm radio wavelengths. 
The bolometric luminosity has contributions from the accretion of infalling material and nuclear 
burning. Up to date only a few systematic searches for high mass protostellar objects have been 
carried out (Molinari et al. 1996, 1998, 2000; Sridharan et al. 2002). 

We have recently started a multi-wavelength study of a sample of 18 luminous IRAS sources in 
the southern hemisphere thought to be representative of young massive star forming regions (Mar- 
dones, Garay, & Bronfman 2002). The goal is to understand the physical and chemical differences 
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between different stages of early evolution. The objects were taken from the Galaxy-wide survey of 
CS(2^1) emission towards IRAS sources with IR colors typical of compact H II regions (Bronfman, 
Nyman, & May 1996). We selected sources based primarily on the observed CS(2^1) line profiles; 
looking for self-absorbed lines consistent with inward or outward motions (e.g., Mardones 1998), 
and/or with extended line wings, possibly indicating the presence of bipolar outflows. In addition, 
the sources were required to have IRAS 100/xm fluxes greater than 10^ Jy and to be in the southern 
hemisphere {5 < —20°). The luminosity of the IRAS sources, computed using the IRAS energy 
distribution and the distances derived by Bronfman (2002) are in the range 2 x 10^ — 4 x 10^ Lq, 
implying that they contain at least an embedded massive star. 

Most of the selected objects are expected to be associated with ultra compact (UC) H II regions 
which are thought to be manifestations of newly formed massive stars that are still embedded in 
their natal molecular cloud. This expectation is confirmed by the radio continuum observations 
of Walsh et al. (1998) which show that 9 of the 12 sources in both samples have detectable radio 
continuum emission (above a 3a limit of 6 mJy/beam at 8.64 GHz with angular resolution of ~1.5"). 
The objects that were not detected at radio wavelengths are suitable candidates for massive stars 
in very early stages of evolution in which dense material is still falling toward a massive OB-type 
protostar. In this accretion phase, the high-mass accretion rate of the infalling material quenches 
the development of an UC H II region (Yorke 1984; Walmsley 1995), and the free- free emission from 
the ionized material is undetectable at centimeter wavelengths. The mass accretion rate might also 
be large enough that the ram pressure of the infalling gas could provide the force to prevent the 
expansion of an H II region. We note, however, that due to the limited sensitivity of the Walsh et 
al. (1998) survey, low-density H II regions with emission measures smaller than 4 x 10^ pc cm~^ 
were not detectable. Hence, the lack of detection at the above limit does not rule out the presence 
of an optically thin compact H II region. 

In this paper we report sensitive ATCA radio continuum observations toward two sources in our 
sample, IRAS 15596—5301 and 16272—4837, without previously detected radio continuum emission 
to place stringent limits in their radio flux density. The la sensitivity level of 70 /xJy at 4.8 GHz 
achieved in our observations is thirty times smaller than in previous studies and is sufficient to detect 
the emission measure corresponding to any ionizing OB star within the Galaxy. The main goal was 
to establish whether or not these objects correspond to very young massive objects; that is massive 
protostars which have not yet ionized significant amounts of the surrounding gas. Also reported 
in this paper are millimeter continuum and molecular line observations of IRAS 15596—5301 and 
16272—4837 made with the Swedish-ESO submillimetre telescope. The latter observations are part 
of a molecular line survey toward several high-mass star forming regions made in order to determine 
their physical characteristics and investigate possible chemical differences. 
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2. 



Observations 



The observations were made using the Austraha Telescope Compact Array ATCA in Aus- 
traha, and the 15 m Swedish-ESO Submillimetre Telescope (SEST) located on La Silla, Chile. 



The ATCA radio continuum observations were made in 2000 May, 23, using the 6D configura- 
tion, which utilizes all six antennas and covers east-west baselines from 77 m to 5.9 km. Observations 

were made simultaneously at the frcqiicncies of 4.800 and 4.928 GHz, each spanning a bandwidth 
of 128 MHz. At these frequencies the FWHM primary beam of ATCA is 10'. Total integration 
time for each source was 270 minutes, obtained from 15-minute scans taken over a wide range of 
hour angles to provide good (u,v) plane coverage. The calibrator PKS 1600-48 was observed before 
and after every on-source scan in order to correct the amplitude and phase of the interferometer 
data for atmospheric and instrumental effects. The flux density was calibrated by observing PKS 
1934-638 (3C84) for which values of 5.83 Jy at 4.8 GHz and 5.65 Jy at 4.9 GHz were adopted. 
Standard calibration and data reduction were performed using MIRIAD (Sault, Teuben, &: Wright 
1995). Maps were made by Fourier transformation of the uniformly weighted interferometer data 
using the AIPS task MX. The noise level in each image is 70 /nJy, close to the expected theoretical 
limit of 50 /xJy. The synthesized (FWHM) beams were 2.30" x 1.68" for IRAS 15596-5301 and 
2.46" X 1.64" for IRAS 16272-4837. 



The 1.2 mm continuum observations were made in 2001, September using the SIMBA 37- 
channel bolometer array. The HPBW of a single element is 24" and the separation between elements 
on the sky is 44". We made the observations in the fast mapping mode, using a scan speed of 80" 
s~^. Our observing blocks consisted of 50 scan lines in azimuth of length 800" and separated 
in elevation by 8", giving a map size in elevation of 400". This block required '^15 minutes of 
observing time. Two blocks were observed for each source. The data were reduced in a standard 
fashion, applying sky-opacity correction, calibration using a counts-to-flux conversion factor derived 
from maps of bright planets, baseline subtraction, and rejection of correlated sky-noise using the 
software package MOPSI. The rms noise level in the final maps is 0.1 Jy with a pixel binning of 8". 
Uncertainties in the pointing accuracy are estimated at la = 5". 



2.1. 



ATCA 



2.2. 



SEST 



2.2.1. Millimeter continuum 



^Tho Australia Telescope Compact Array is funded by the Commonwealth of Australia for operation as a National 
Facility managed by CSIRO 
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2.2.2. Molecular lines 

The molecular line observations were carried out in two epochs. During June of 1999 we used 
as backend the high resolution acousto-optical spectrometers which provided a channel separation 
of 43 kHz and a total bandwidth of 43 MHz. We observed, at the position of the IRAS sources, 
molecular line emission in ten transitions: CS(2^1), CS(5->4), CS(7^6), C3^S(2^1), SiO(2^1), 
SO(65 ^ 54), HCO+(1^0), HCO+(3^2), Hi3C0+(l ^0) and 0^^0(2^!), with integration times 
on source ranging from 3 to 6 minutes. System temperatures were typically ~ 220 K at 3 mm, 
^ 350 K at 1 mm, and ~ 700 K at 0.8 mm. In addition we mapped, with angular spacings of 
30", the CS(2^1) and CS(5^4) emission within regions of 2.5' and the SiO(2^1) emission within 
regions of 2.5' x 1.5' centered on the IRAS positions. These observations were made in order to 
investigate the kinematics of the molecular gas. 

During March of 2000 we made a survey of molecular line emission using as backend the low 
resolution acousto-optical spectrometers providing a channel separation of 700 kHz and a total 
bandwidth of 1.0 GHz. The 100/150 GHz and 115/230 GHz pairs of SIS receivers were used to 
simultaneously observe lines in, respectively, the 3/2 mm and 3/1 mm wavelength ranges. We 
observed 39 spectral windows, 1.0 GHz wide each, within the 85 to 250 GHz frequency range. 
System temperatures were typically ^ 200 K at 2-mm and 3-mm, and ~ 400 K at 1-mm. The 
integration times on source were 5 minutes, resulting in antenna temperature rms noises of typically 
20 mK. A detailed description of these observations are given in Garay et al. (2002). All the 
observations were performed in the dual beam-switching mode, with a beam separation of 11'47" 
in azimuth. 



3. Results 

3.1. Radio continuum emission. 

Figure 1 shows a map of the 6 cm radio continuum emission from IRAS 15596—5301. Three 
compact radio sources (labeled A, B, and C) were detected within a region of 30" in diameter. 
Their positions, flux densities, and sizes are given in Table 1. Their peak flux densities are in the 
range 1.1 to 2.8 mJy/beam, close to the sensitivity limit in the survey of Walsh et al. (1998). The 
less compact source (object A) appears to have an irregular shell-like structure. 

No radio continuum emission was detected toward 16272-4837 within a region of 3' in radius 
centered on the IRAS position. We set a new 3a upper limit for the 6 cm flux density of 16272-4837 
of 0.2 mJy, about thirty times smaller than that previously reported by Walsh et al. (1998). We 
detected, however, emission from a region located about 4' northwest of the array phase center, 
with a total flux density (corrected by primary beam) of 122 mJy. The morphology of the emission 
is complex, showing an elongated bipolar-like structure with a peak located at 0:2000 = 16'^30"^35?0, 
^2000 = -48°40'47.6". 
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3.2. Millimeter continuum emission. 

Figure 2 presents maps of the 1.2-mm continuum emission observed with SEST, revealing the 
presence of strong emission towards both IRAS sources. The 1.2-mm emission associated with 
IRAS 15596-5301 has a peak position at q;2ooo = 16''03™32!9, (52ooo = -53°09'20", a flux density 
of 5.8 ± 0.3 Jy and a peak flux density of 2.0 Jy/bcam. It shows an elongated morphology, with 
deconvolved major and minor axes of 42.1" and 24.6" at P. A. 46°. The 1.2-mm emission detected 
towards IRAS 16272-4837 has a peak position at 02000 = 16'^30™58!7, (^2000 = -48°43'55", a flux 
density of 13.8 ± 0.9 Jy and a peak flux density of 4.6 Jy/beam. It has an elongated morphology, 
with deconvolved major and minor axes of 41.0" and 24.7" at P.A. 142°. 

3.3. Molecular line emission. 

The molecular line survey indicate that IRAS 15596—5301 and 16272—4837 have similar chem- 
ical characteristics. This is illustrated in Figure 3 which shows rotational diagrams of the emission 
from selected species (CH3OH, CH3CN, CH3CCH, and HCCCN) with optically thin hues and at 
least five observed transitions. The average rotational temperature derived from these diagrams is 
27 K for both sources. Overall wc detected emission from more than 20 species toward both objects. 
Table 2 summarizes the derived properties of the molecular emission in selected species. Column 
densities, rotational temperatures, and abundances relative to CO are given in columns 2-4. The 
column density of optically thick species, such as CO and CS, were determined as follows. The 
optical depths of a pair of an optically thick line and an optically thin line of an isotopic species are 
determined, for a given abundance ratio, from the ratio of their observed brightness temperatures. 
From these opacities and assuming an excitation temperature of 30 K we then computed the total 
column density assuming that energy levels are populated according to local thermodynamic equi- 
librium (see Bourke et al. 1997 for a description of the method). For the derivation of CO column 
densities we used the optically thick ^^C0(1 —>■ 0) and optically thin C^^0(1 0) pair of lines, 
whereas for CS we used the CS(2 1) and C^^S(2 1) pair of lines. For species with optically 
thick lines it is also possible to determine the filling factor of the emission within the observed 
beam. Filling factors of ~0.5 for ^^CO and ~0.13 for CS were estimated for both sources. 

Figure 4 shows maps of the CS(5 4) line emission from both sources. The morphology of the 
dense gas from IRAS 15596—5301 is roughly circular, with a mean angular diameter of 37". This 
implies a cloud radius of 0.41 pc at a distance of 4.6 kpc. The spatial distribution of the dense gas 
from IRAS 16272—4837 is elongated, with major and minor axes of 61" and 42", respectively. This 
implies radii of 0.50 and 0.35 pc at a distance of 3.4 kpc. Prom the observed CS size and derived 
CS column density we estimate that the CS density is 4 x 10~^ cm~^ in IRAS 15596—5301 and 
2 X 10~^ cm~^ in 16272—4837. Further, assuming a [CS/H2] abundance ratio of 1 x 10^^ (cf. van 
der Tak et al. 2000), we find that the molecular hydrogen density and the mass of molecular gas 
are 4 x 10^ cm-^ and 6.6 x 10^ Mq for IRAS 15596-5301, and 2 x 10^ cm'^ and 2.6 x 10^ Mq for 
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IRAS 16272—4837. Thus, the molecular observations show that both sources are associated with 
dense, massive molecular cores. The mass of the cores can, alternatively, be estimated assuming 
that they are in virial equilibrium (MacLaren et al. 1988). From the observed size and average line 
width in the CS(5-4) line we derive virial masses of 2.3 x 10^ M© for IRAS 15596-5301 {R = 0.41 
pc; Ai;=5.17 km s-^) and 2.6 x 10^ Mq for IRAS 16272-4837 {R = 0.42 pc; A?;=5.44 km s^^). 

Figure 5 shows the spectra observed, with the high resolution spectrometer, toward the peak 
position of the IRAS 16272-4837 core. The CS(2^1), HCO+(1^0) and HCO+(3^2) spectra show 
double-peaked line profiles, with a bright blue-shifted peak at the velocity of ~ —47.1 km s~^ and 
a weaker red-shifted peak at the velocity of ~ —44.6 km s"-*^. On the other hand, the profiles of the 
C^^S(2^1) and H^^C0"*"(1— >0) lines show a symmetric single component with a peak line center 
velocity of '-^ —46.0 km s~^. These spectroscopic signatures suggest that the bulk of the molecular 
gas toward IRAS 16272—4837 is undergoing large-scale inward motions. Infalling motions traced 
by optically thick molecular lines are expected to produce hne profiles showing blue asymmetry, 
whereas optically thin lines are expected to exhibit symmetrical profiles (Mardones et al. 1997). 

In some species, particularly SiO and SO, the spectra show the presence of strong wing emission 
indicative of outflow activity. In the SiO (2 — > 1) line the wing emission extends up to velocities of 
—19.6 km s"^ and 19.5 km s~^ relative to the ambient cloud velocity of —46.2 km s~^. Figure 6 
shows contour maps of the velocity integrated SiO emission in the range of flow velocities from 
— 19.5 to —3.8 km s~^ (continuous lines; blueshifted gas) and from 4.2 to 15.9 km s^^ (dotted lines; 
redshifted gas). The flow velocity is defined as Vfiow = vlsr ~ Vq, where Vo is the velocity of the 
emission from quiescent ambient gas. The spatial distribution of the wing emission clearly shows a 
bipolar structure with redshifted and blueshifted lobes extending in opposite directions of a 22/xm 
source detected by MSX (see later discussion). The bipolar structure shows a moderate degree of 
collimation, with the blueshifted emission seen mainly toward the west and the redshifted emission 
seen mainly toward the east, suggesting that the underlying wind is a wide angle wind. 

3.4. Mid-infrared emission. 

We have used the database from the Midcourse Space Experiment (MSX; Price 1995) survey 
of the Galactic plane (Egan et al. 1998) to investigate the mid-infrared emission, in the four MSX 
broad bands (A:6.8-10.8 yum, 0:11.1-13.2 fim, D:13.5-15.9 fim, and E:18.2-25.1 /xm), from both 
IRAS sources. Figure 7 presents grey scale images of the emission in the A, C, and E bands of 
a region of about 4' toward IRAS 15596—5301. These images show the presence of two sources 
within the central 1' region: a compact object, brighter in the E band image, and an extended 
object, brighter in the C band, with a peak located about 30" southwest of the former. Both 
sources are seen in emission in the four MSX bands. The three radio continuum sources detected 
within this region (delineated by the 0.5 mJy contour level in the upper panel of Figure 7), as well 
as the peak position of IRAS 15596—5301 (marked by the cross), lie projected within the face the 
compact object. Further, the extent of the 1.2-mm dust continuum emission (delineated by the 0.5 
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Jy contour level in the lower panel of Figure 7) is similar to that of the E-band infrared emission. 
These results show that the compact MSX object is intimately associated with the 15596—5301 
massive molecular core. We measured flux densities of 5.4, 4.1, and 24.3 Jy in the C, D, and E 
bands, respectively. The extended object shows strong emission in the A and C bands, which are 
sensitive to emission from polycyclic aromatic hydrocarbon (PAH) features (7.7 and 8.3 in A 
band; 11.3 and 12.6 jum in C band). Since no radio emission is detected toward this MSX source, 
we suggest that it might correspond to a photodissociated region excited by a non-ionizing star. 

IRAS 16272—4837 is associated with an infrared dark cloud seen in absorption against bright 
mid-infrared emission in the A, C, and D bands. This is shown in the top two panels of Figure 8 
which presents grey scale images of the emission in the A, C, and E bands toward IRAS 16272—4837. 
There is a good agreement between the molecular, millimeter and mid-infrared morphologies, as 
illustrated by the contour of the 1.2-mm continuum emission shown in Figure 8. From the molecular 
hydrogen column density derived from the line observations (of 3.4 x 10^^ cm~^), assuming a gas-to- 
dust mass ratio of 110 (Draine & Lee 1984) and a dust opacity at 8/Ltm of 10^ cm^ gr"^ (Ossenkopf 
& Henning 1994), we estimate that the optical depth at 8/xm is ~ 12, indicating that this cloud is 
optically thick at MID-IR. Further, the MSX observations in the E-band show the presence, near 
the center of the infrared dark cloud, of an object seen in emission located at 02000 = 16'^30™58!18, 
(^2000 = -48°43'47.9" (see bottom panel of Figure 8). We finally note that the MSX object seen 
near the top right corner of Figure 8, which is clearly extended in the A-band image, is associated 
with the IRAS object 16269-4834 and coincident with the radio source detected ^4' northwest of 
16272-4837. 



4. Discussion 

4.1. Spectral energy distribution 

Figure 9 shows the spectral energy distribution (SED) of IRAS 15596—5301 and 16272—4837 
from 12 fim to 1.2 mm, which is mainly due to thermal dust emission. We fitted the SED with 
modified blackbody functions of the form Bi,{Tcl) [1 ~ exp(— r,^)] Q,s , where is the dust optical 
depth, By{Td) is the Planck function at the dust temperature T^, and fi^ is the solid angle subtended 
by the dust emitting region. The opacity was assumed to vary with frequency as z/^, i.e. = 
{u/Ufj)^ , where Vo is the frequency at which the optical depth is unity. Due to the limited number 
of spectral points we have set the value of /3 equal to 2.0, consistent with tabulated opacities 
(Ossenkopf k, Henning 1994) and derived values for high mass star forming regions (Molinari et al. 
2000). A single temperature model produced poor fits, underestimating the emission observed at 
wavelengths smaller than 25/xm, and therefore we used a model with two temperature components. 

From the fits (long-dashed lines) we derive that the colder dust component (short-dashed lines) 
toward 15596—5301 and 16272—4837 have, respectively, temperatures of 27 and 25 K, angular sizes 
(assuming a Gaussian flux distribution) of 30"(FWHM), and wavelengths at which the opacity is 
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unity of ~90 fim and 140 /xm. The temperature of the hot dust component is 100 K for 15596—5301 
and 115 K for 16272—4837. The thermal dust emission at 1.2 mm is therefore optically thin 
(7"i.2mm 5 X 10~^), and thus the observed flux density at 1.2 mm allows to obtain an additional 
mass estimate of the dense cores. In general, for an isothermal dust source the total gas mass, Mg, 
is given in terms of the observed flux density, S^, at an optically thin frequency, v, by (e.g. Chini, 
Krugel, & Wargau 1987) 

M = " 

where is the mass absorption coefficient of dust, R^g is the dust-to-gas mass ratio (assuming 
10% He), and By{Td) is the Planck function at the dust temperature T^. The main source of 
uncertainty in the conversion of the observed flux density into gas mass is the RdgHu factor, or 
total mass opacity, which is a poorly known quantity (e.g. Gordon 1995). Using a dust opacity at 
1.2 mm of 1 cm^ g~^, as computed by Ossenkopf &; Henning (1994) for dense and cold protostellar 
cores, Rdg = 0.01, the fltted dust temperatures, and the observed flux densities, we derive masses 
of 1.4 X 10^ Mq for IRAS 15596-5301 and 2.0 x 10^ M© for IRAS 16272-4837. These masses 
derived from the dust emission are in good agreement with those derived from the molecular line 
intensities and from the virial assumption. 



4.2. Evolutionary stages 

4.2.1. IRAS 15596-5301 (GS29.40-0.46) 

The radio continuum observations toward IRAS 15596—5301 indicate the presence, in a region 
of ~0.3 pc in radius, of three distinct H II regions with diameters ranging from 0.06 to 0.2 pc. The 
multiple structure of the ionized gas is typical of galactic H II regions, and is most likely due to the 
presence of a cluster of exciting stars. If components A, B, and C are excited by individual ZAMS 
stars, the rate of UV photons needed to ionize them (see Table 3) imply exciting stars with spectral 
types of BO, BO. 5, and Bl, respectively. The total luminosity emitted by this cluster of B stars, 
as inferred from the radio observations, is 4.1x10^ Lq. On the other hand, the total far-infrared 
luminosity computed using the IRAS fluxes (see Casoli et al. 1986) is 6.5 x lO'^ Lq (assuming a 
distance of 4.6 kpc; Bronfman 2002). The difference between the radio derived luminosity and the 
IRAS luminosity could be explained by the presence of dust within the H II regions. Garay et al. 
(1993) found that the fraction of Lyman continuum photons absorbed by dust within H II regions 
is typically 55%. Alternatively, it could be explained by the presence, in addition to the B stars, of 
several less massive stars that will contribute to the FIR luminosity but that are not hot enough 
to contribute to ionization. It is not easy, however, to disentangle which of these effects play the 
predominant role. 

The compact H II regions are found projected toward the peak of the CS(5 4) emission map 
(see Figure 4), suggesting that they are deeply embedded within the dense molecular core. From 
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the observed sizes, and assuming a sound speed in the ionized gas of 11.4 km s~ , we estimate 
that the H II regions have dynamical ages between 3x10^ to 8x10^ yrs. If these correspond to the 
actual ages of the compact H II regions, then we should conclude that they are very young objects. 
The dynamical time-scales, however, may not provide a realistic estimate of the actual age of H II 
regions. The large number of UC H II regions and their short dynamical ages poses the well known 
problem that the rate of massive star formation appears to be much greater than other indicators 
suggest (Wood & Churchwell 1989, Churchwell 1990). Due to the high density of the molecular gas 
in which they are embedded, we suggest instead that the H II regions within the IRAS 15596—5301 
massive core might be in pressure equilibrium with the surrounding dense ambient medium, and 
are currently stalled at their equilibrium radius, (e.g. De Pree, Rodriguez, & Goss 1995). The 
molecular density of the ambient gas needed to stall an H II region at radius, Rf, is given by (e.g. 
Garay &; Lizano 1999) 




where is the rate of ionizing UV photons emitted by the exciting star, Tg is the electron 
temperature of the ionized gas, and Tq is the temperature of the ambient gas. Using the observed 
radius and the derived ionizing rate of UV photons of the H II regions within IRAS 15596—5301, we 
find that molecular densities of ~ 5 x 10^ cm~^ are needed for them to be pressure confined by the 
dense environment. These densities are similar to those derived from the molecular observations. 
The time needed for the H II regions to achieve pressure equilibrium are between 1.0 x 10^ yrs to 
2.5 X 10^ yrs, implying that massive star formation started within this core more than 2.5 x 10^ yrs 
ago. We conclude that the dense massive core is in an advanced stage of early evolution, in which 
multiple OB star formation have already taken place near its central region. 



4.2.2. IRAS 16272-4837 (G335. 58-0.28) 

The total far-infrared luminosity of IRAS 16272—4837 computed using the IRAS fluxes (see 
CasoU et al. 1986) is ~ 2.4 x 10^ Lq (assuming a distance of 3.4 kpc; Bronfman 2002). The 

luminosity obtained integrating under the fitted curve in Figure 9 is similar to the IRAS luminosity. 
The high luminosity suggests that the IRAS 16272—4837 massive core hosts a young massive 
protostar inside; and the lack of radio emission suggests that it is still undergoing an intense 
accretion phase. Models of massive envelopes accreting onto a young massive central B type star 
(e.g. Osorio et al. 2000) require accretion rates in the envelopes M > 5 x 1O~^M0 yr~^ in 
order to fit the observed SEDs. The high-mass accretion rate of the infalling material quenches 
the development of an UC H II region (Yorke 1984), and the free- free emission from the ionized 
material is undetectable at centimeter wavelengths. 

This hypothesis is supported by the characteristics of the line profiles observed toward IRAS 
16272—4837, which suggest that the molecular gas is undergoing infalling motions. Prom the spectra 
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of the optically thick HC0"'"(1 —>■ 0) line we measure a velocity difference between the red and blue 
peaks of 2.7 km s~^, and brightness temperature of the blue peak, red peak, and dip of 5.4, 3.6, 
and 2.9 K, respectively. From the spectra of the optically thin H^'^C0+(1 0) line we measure a 
FWHM line width of 3.18 km s~^. From these values, using the simple model of contracting clouds 
of Myers et al. (1996), we derive a characteristic inward speed of 0.5 km s~^. We note that this 
value is considerably smaller than the free-fall velocity expected for a cloud with a total mass of 
~ 3 X W^Mq at its outer envelope radius of 0.4 pc, suggesting that the collapse is not dynamical. 
Using the derived values of the infall speed, molecular density, and core size, we obtain a mass infall 
rate, Mj„, of 1 x IQ-^Mq yi'^, large enough to prevent the development of an UC H II region. 
The high value of the mass infall rate rises the question as to which is the fraction of the total 
luminosity due to accretion. The accretion luminosity, Lace, is 

_ GfMin Mp 
^acc — T-, ) 

Up 

where / is the fraction of the large scale mass infall rate that goes into accretion onto the protostar, 
Mp is the mass of the protostar, and Rp the radius where the shock occurs. None of these three 
parameters are known for IRAS 16272-4837. Assuming Mp ~ IOMq, ~ 3 x 10^^ cm, and 
/ ~ 0.05 (e.g, Norberg & Maeder 2000) we obtain Lace ~ 3.6 x 10^ Lq, about 15% of the total 
luminosity. We emphasize that this value of the accretion luminosity corresponds only to a rough 
estimate, particularly because the value of / is highly uncertain. 

Additional evidence for IRAS 16272—4837 to be in a collapsing stage is provided by the presence 
of bipolar outflowing gas. phenomenon which is thought to be closely related to accretion processes. 
The 22fim object lies at the center of symmetry of the SiO outflow, suggesting it is intimately 
associated with the energy source of the outflow. The early evolutionary stage of this region is 
also sustained by the presence of 6.67 GHz methanol masers, which are thought to be signposts of 
young regions of massive star formation (Walsh et al. 1997, 1998). Notice that the 22/xm source 
is elongated, with the 6.67 GHz methanol masers being aligned along its major axis. Finally, we 
mention that the high value of the mass to luminosity ratio of IRAS 16272—4837, M/L = 0.083, 
about 4 times higher than that of IRAS 15596—5301, is another indicator of its youth, as argued 
by Sridharan et al (2002). 



5. Summary 

We undertook sensitive radio continuum observations at 4.8 GHz, using ATCA, and 1.2 mil- 
limeter continuum and molecular line observations, using SEST, toward IRAS 15596—5301 and 
16272—4837, two luminous objects (£ > 2 x W^Lq) thought to represent massive star forming 
regions in early stages of evolution. The main results and conclusions are summarized as follows: 

1. The 4.8 GHz radio continuum observations toward IRAS 15596—5301 show the presence of 
three distinct compact sources with diameters of 0.06 to 0.2 pc, all located within a region of 30" in 
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diameter. If they are regions of ionized gas, they have electron densities of 8 x 10^ — 2 x 10^ cm~^ 
and are excited by early B type stars. The 1.2-mm continuum observations show that the dust 
emission arises from an elongated region with major and minor axes of 42.1"and 24.6", respectively. 
The observed 1.2-mm flux density, of 5.8 Jy, implies a mass of 1.4 x 10^ Mq. The line observations 
indicate a region of molecular gas with a radius of ~ 0.4 pc, a molecular hydrogen density of 4 x 10^ 
cm~^, and a rotational temperature of ~ 27 K. We conclude that IRAS 15596—5301 corresponds 
to a dense massive (~ 3 x 10^ -^0) molecular core, with a radius of about 0.4 pc, that is in the 
compact H II region phase of evolution, already hosting a cluster of massive B type stars. Further, 
we suggest that the regions of ionized gas excited by these stars are in pressure equilibrium with 
the ambient molecular gas. 

2. No radio continuum emission at 4.8 GHz was detected from IRAS 16272—4837 up to a 3a 
limit of 0.2 mJy. The 1.2-mm continuum observations show strong dust emission arising from a 
region with an elongated morphology, with major and minor axes of 41.0" and 24.7". The observed 
flux density at 1.2-mm, of 13.8 it 0.9 Jy, implies a mass of 2.0 x 10^ Mq. The line observations 
indicate a molecular gas region with a size of ~ 0.4 pc, a molecular hydrogen density of 2 x 10^ 
cm~^, and a rotational temperature of ~ 27 K. We suggest that this object corresponds to a dense 
massive core in a very early evolutionary stage, distinguished by being luminous but not associated 
with an UC H II region (eg. Cesaroni et al. 1994; Hunter et al. 1998; Molinari et al. 1998). The 
core is being heated by a recently formed massive star embedded at its center which is still accreting 
at high rates, the main heating agent being the accretion luminosity, and have not yet produced 
a detectable ultra compact H II region. As such, this will become a key object for the study at 
submillimeter wavelengths with high angular resolution of the earliest stages of the formation of 
massive stars, and their related accretion and outflow processes. 

G.G. and D.M. gratefully acknowledge support from the Chilean Ccntro de Astrofisica FON- 
DAP No. 15010003. M.G.B. acknowledges support from the Australian Research Council. 
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Table 1. OBSERVED PARAMETERS OF RADIO CONTINUUM EMISSION 



H II region Position Flux density Angular size 

a (2000) (5(2000) (mJy) ("x") 



A 16''03"'31f01 -53° 09' 32:'9 43.4 9.7x8.0 

B 16 03 31.80 -53 09 21.5 9.0 3.0x2.5 

C 16 03 31.91 -53 09 31.1 3.0 6.1x2.6 



Table 2. MOLECULAR GAS PARAMETERS 



Molecule IRAS 15596 IRAS 16272-4837 

N [X/CO] N [X/CO] 

(K) (cm-2) (K) (cm-2) 



CO 


30'' 


8.1 


x 


1019 




1 




30^ 


2.1 


X 


1019 




1 




i^CO 


30^ 


1.5 


X 


IQis 


1.9 


X 10" 


-2 


30^ 


3.8 


X 


IQi^ 


1.8 


X 10" 


-2 


C^^o 


30^ 


1.1 


X 


IQi^ 


1.4 


X 10" 


-3 


SO'^ 


2.7 


X 


1016 


1.3 


X 10" 


-3 


cs 


30^ 


9.8 


X 


1015 


1.2 


X 10^ 


-4 


30^ 


4.2 


X 


1015 


2.0 


X 10" 


-4 


C34S 


30^ 


4.4 


X 


10^4 


5.4 


X 10- 


-6 


30^ 


1.9 


X 


1014 


9.0 


X 10" 


-6 


H2C0 
















27.0 


9.4 


X 


1013 


4.5 


X 10- 


-6 


CH30H 


19.1 


2.7 


X 


10^4 


3.3 


X 10" 


-6 


18.0 


5.5 


X 


1014 


2.6 


X 10" 


-5 


CH3CCH 


26.3 


2.1 


X 


1014 


2.6 


X 10" 


-6 


31.9 


2.3 


X 


1014 


1.1 


X 10" 


-5 


CH3CN 


33.7 


9.0 


X 


10^2 


1.1 


X 10" 


-7 


29.0 


1.9 


X 


1013 


9.0 


X 10" 


-7 


CH3CHO 


25.0 


7.2 


X 


10^2 


8.9 


X 10" 


-8 


21.0 


1.3 


X 


1013 


6.2 


X 10" 


-7 


SiO 


10.4 


2.9 


X 


10^2 


3.6 


X 10- 


-8 


9.9 


1.0 


X 


1013 


4.8 


X 10- 


-7 


SO 


23.0 


6.3 


X 


1013 


7.8 


X 10- 


-7 


11.7 


8.5 


X 


1013 


4.0 


X 10- 


-6 


HCCCN 


28.2 


4.1 


X 


10^3 


5.1 


X 10" 


-7 


27.3 


5.1 


X 


1013 


2.4 


X 10- 


-6 


OCS 


37.5 


3.9 


X 


1013 


4.8 


X 10" 


-7 


28.5 


1.3 


X 


1014 


6.2 


X 10- 


-6 


SO2 


27.6 


1.9 


X 


1013 


2.3 


X 10" 


-7 


25.9 


1.6 


X 


1013 


7.6 


X 10- 


-7 



^Assumed excitation temperature. 



Table 3. DERIVED PARAMETERS OF H II REGIONS 



H II region 


Diameter 


EM 


Tie 


N„ 


S.T. 




(pc) 


(pc cm~^) 


(cm^^) 


(s-i) 




A 


0.196 


2.4 X 10^ 


9.2 X 102 


8.3 X lO''^ 


BO 


B 


0.061 


5.2 X 10^ 


2.4 X 10^ 


1.7 X 10^6 


B0.5 


C 


0.089 


8.0 X 10"^ 


7.8 X 10^ 


5.7 X 10^5 


Bl 
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Fig. 1. — ATCA map of the 6-cm radio continuum emission from IRAS 15596—5301. The angular 
resolution is 2.30" x 1.68" (FWHM beam is shown in the lower left corner). Contour levels are -3, 
3, 6, 9, 12, 18, 24, 30, and 36 times 70 /xJy beam"^ (= la noise level). 
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Fig. 2. — Maps of 1.2-mm continuum emission observed with SEST. The angular resolution is 

24"(FWHM beam is shown in the lower left corner). Top: IRAS 15596—5301. Contour levels are 
-1, 1, 2, 3, 4, 5, and 6 times 0.3 Jy beam'^ Bottom: IRAS 16272-4837. Contour levels are -1, 1, 
2, 3, 4, 5, 6, 7, 8, 9 and 10 times 0.5 Jy beam^^ 
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Fig. 3. — Rotational diagrams of the emission from selected species observed toward the IRAS 
15596-5301 and 16272-4837 massive dense cores. From top to bottom: CH3OH, CH3CN, 
CH3CCH, and HCCCN. The lines correspond to least squares linear fits to the observed data. 
The derived rotational temperature is given in the upper right corner. 
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- 15596-5301 CS (5-4) - 




Right Ascension (J2000) 

Fig. 4. — Maps of CS(5^4) line emission observed with SEST. The angular resolution is 
22"(FWHM beam is shown in the lower right corner). Top: IRAS 15596—5301. Contour lev- 
els are 20, 30, 40, 50, 60, 70, 80, and 90 % of the peak velocity integrated antenna temperature 
of 11.4 K km s~^. The triangles mark the positions of the compact H II regions. Bottom: IRAS 
16272-4837. Contour levels arc 20, 30, 40, 50. 60, 70, 80, and 90 % of the peak velocity integrated 
antenna temperature of 7.2 K km s^^. The star marks the peak position of the 22/xm MSX source. 
In both panels the cross mark the position of the IRAS source and the squares the position of 6.67 
GHz methanol masers (Walsh et al. 1998). 
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Fig. 5. — Spectra of the emission in several molecular lines observed toward the peak position of 
the IRAS 16272—4837 massive core. Transitions are given in the upper right corner of the spectra. 
The vertical dotted line indicates the systemic velocity of the ambient gas of —46.2 km s~^. 
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Fig. 6. — Contour maps of velocity integrated SiO (2 1) line wing emission from IRAS 
16272—4837. Continuous lines represent blue-shifted emission {vfio-w ■ —19.5 to —3.8 km s~^) 
and dashed lines red-shifted emission {vfiow • 4.2 to 15.9 km s~^). The star marks the position of 
the 22/xm MSX source. Contour levels are at 0.60, 0.72, 0.84, 0.96, 1.08 and 1.20 K km s^^. 
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Fig. 7. — MSX images (gray-scale) towards IRAS 15596—5301. The cross marks the peak position 
of the IRAS source. Top: A-band image. The contour dchneates the radio continuum emission at 
the level of 0.5 mJy. Middle: C-band image. The squares and triangles mark the position of 6.67 
GHz methanol and hydroxyl masers, respectively. Bottom: E-band image. The dotted contour 
indicates the 1.2-mm dust emission at the level of 0.5 Jy. 




Fig. 8. — MSX images (gray-scale) towards IRAS 16272—4837. The cross marks the peak position 
of the IRAS source and the contour dehneates the 1.2-mm continuum emission at the level of 0.5 
Jy. Top: A-band image. The light areas in the gray-scale correspond to regions of lower emission 
than the bright Galactic plane emission, and indicates the infrared extinction cloud associated with 
IRAS 16272—4837. Middle: C-band image. The squares and triangles mark the position of 6.67 
GHz methanol and hydroxyl masers, respectively. Bottom: E-band image. 
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Fig. 9. — Spectral energy distributions. Stars mark IRAS fluxes, squares MSX fluxes, and the circle 
the SIMBA flux. The long-dashed curve is a fit to the spectrum using two modified blackbody 
functions of the form B^{T^) [l — cy:p{—{v/vo)^)], with different temperatures. The short-dashed 
line indicates the fit for the colder temperature component (fit parameters indicated on the upper 
left). Top: IRAS 15596-5301. Bottom: IRAS 16272-4837 (12/xm flux is an upper limit). 



